Draft version February 5, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



THE INFRARED PROPERTIES OF SUBMILLIMETER GALAXIES: CLUES FROM ULTRA-DEEP 70 /im 

IMAGING 

MlNH T. HUYNH 

Spitzer Science Center, MC220-6, California Institute of Technology, Pasadena CA 91125, USA. mhuynh@ipac.caltech.edu 

Alexandra Pope 

Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, V6T 1Z1, Canada 



O 
O 

(N 

O 

<D 

Q 

o 

(N 

(N 
> 
O 

in 

(N 
O 

Oh: 

6 



David T. Frayer 

Spitzer Science Center, MC220-6, California Institute of Technology, Pasadena CA 91125, USA 

Douglas Scott 

Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, V6T 1Z1, Canada 

Draft version February 5, 2008 

ABSTRACT 

We present 70 /im properties of submillimeter galaxies (SMGs) in the Great Observatories Origins 
Deep Survey (GOODS) North field. Out of thirty submillimeter galaxies (Ssso > 2 mJy) in the central 
GOODS-N region, we find two with secure 70 /im detections. These are the first 70 /Ltm detections of 
SMGs. One of the matched SMGs is at z ~ 0.5 and has S-jq/ S^q and <SVo/ '£24 ratios consistent with a 
cool galaxy. The second SMG (z = 1.2) has infrared-submm colors which indicate it is more actively 
forming stars. We examine the average 70 /xm properties of the SMGs by performing a stacking 
analysis, which also allows us to estimate that Sgso > 2 mJy SMGs contribute 9 ± 3% of the 70 /im 
background light. The Ssso/SVo colors of the SMG population as a whole is best fit by cool galaxies, 
and because of the redshifting effects these constraints are mainly on the lower z sub-sample. We fit 
Spectral Energy Distributions (SEDs) to the far-infrared data points of the two detected SMGs and 
the average low redshift SMG (z mc dian = 1.4). We find that the average low- z SMG has a cooler dust 
temperature than local ultraluminous infrared galaxies (ULIRGs) of similar luminosity and an SED 
which is best fit by scaled up versions of normal spiral galaxies. The average low-z SMG is found to 
have a typical dust temperature T= 21-33 K and infrared luminosity is-iooOjum = 8.0 x 10 11 L Q . We 
estimate the AGN contribution to the total infrared luminosity of low-z SMGs is less than 23%. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: starburst — submillimeter 



1. INTRODUCTION 

Deep submillimeter surveys provide a probe of galax- 
ies which is almost independent of luminosity for a 
wide redshift range 1 < z < 5, beca use of the negativ e 
K-correction in the infrared (e.g. iBlain et al.l [2002). 
The Sub millimeter Common User Bolometer Array 
(SCUBA; iHolland et all fl999h and Max-Planck Mil- 



limeter Bolometer (MAMBO. iKrevsa et al.lll998| ) have 
now discovered several hundred s ubmillimeter galax- 
ies (SMGs: e.g. iSmail et all ll997L iHughes et al.l ll998L 
Barger et all 119981. iBorvs et alJ 120021 
20021 I Webb et alJ I2003L IGreve et al " ' 



zuuzl ivveuu eu ai.i izuuot ivjieve et ai.i izuu'a i.rupe et ai.i 
20051 ICoppin et al.ll2006l ). Our knowledge of SMGs has 
been hampered by the large beamsize of the submm 
telescopes, which makes it difficult to identify counter- 
parts in the optical and near-infrared. For example, in 
the Great Observatories Origins Deep Survey (GOODS) 
North field there are typically ~ 10 Hubble Space Tele- 
scope (HST) Advanced Camera for Surveys (ACS) opti- 
cal galaxies within a single SCUBA 15 arcsec beam. 

Deep radio surveys have proven to be one of 
the best ways to iden t ify the counterpar t s to SMGs 
fe-g jlvison et alJ I200CI iBarger et ail I2000L iPope et all 
2006). The star formation processes that heat the dust 



responsible for the submm light also produce radio emis- 
sion, as evidenced by the well known far-infrared-radio 
correlation. Taking advantage of this, SMGs have been 
matched to radio sources, which have much better po- 
sitional accuracy, and the optical/near-IR counterparts 
can then be found after the radio identification has been 
made. This is possible because the density of sources in 
the deepest radio images is much less than that of opti- 
cal galaxies, so the probability of a chance coincidence is 
lower. 

The current knowledge of SMGs is biased towards the 
radio identified sources. The SMGs have a median red- 
shift of about 2.2 an d bolometric infrared (IR) lumi- 
nosities Lir > 10 12 L (jChapman et al.ll2004l2005h . The 
optically identified SMGs are faint optically, (1775 >22), 
and red (1775 — K a ~ 2.3 in AB magnitudes), with about 
30% having colors consistent with Extremely Red Ob- 
jects fEROs. lPope et alj|2005fl . Thus SMGs are thought 
to be high redshift, dusty analogs of local ultraluminous 
infrared galaxies (ULIRGs). 

The spectral energy distributions (SEDs) of SMGs 
have traditionally been fit with loc al ULIRGs as tem - 
plates, and in particular Arp220 (e.g. [Barger et alT 2000). 
which has an effective dust temperature of about 42- 
47 K (jKlaas et al.lll997llDunne et al.ll2000f) . Chapman et 
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al. (2005) found a typical dust temperature of 36 ± 7 K 
and a median Lir, = 8.5 x 10 12 L© for their SMG sample. 
This is cooler than local ULIRGs, which have an aver- 
age dust temperature of 43 ± 6 K (based on the Dunne et 
al. 2000 sample). The Far-IR BACKground (FIRBACK) 
study of 170 /im selected galaxies found two ULIRGs be- 
tween 0.5 < z < 1 w hich have SEDs co oler and less lumi- 
nous than Arp220 (jSaiina et al.ll2006h . A 350 /urn study 
of radio-detected SMGs found temperatures of 35 ± 3 K 
(Kovacs et al. 2006). Using the 2 4 /im imaging from the 
Spitzer Legacy Project GOODS, iPope et alj (|2006l ) se- 
curely identify 60% of SMGs in this field in the mid- 
infrared (MIR), and have tentative counterparts for an- 
other 34%. It was found that the observed MIR-submm- 
radio SED of the SMGs peak at longer wavelengths than 
local ULIRGs and are best fit by models with tempera- 
tures of about 30 K (Pope et al. 2006). There is thus an 
emerging picture that SMGs are cooler than previously 
thought. 

Dust temperature affects the inferred IR luminosity, 
and hence star formation rate (SFR), that is derived for 
these galaxies. Studies of SMGs have often assumed dust 
temperatures of 40 K (see lBlain et al.ll2002j ). A drop from 
40 K to 35 K in temperature decreases the IR luminos- 
ity by about a factor of two. Better knowledge of the 
temperature of SMGs is thus crucial for more accurate 
luminosity estimates. 

The advent of the Spitzer Space Telescope makes it 
possible to study the MIR and FIR properties of SMGs 
in detail for the first time. At the median SMG redshift 
of z ~ 2, PAH and silicate features fall into the rest-frame 
wavelength of the 24 /im band. This can make it difficult 
to determine the total infrared luminosities and to fit 
model SEDs well to the FIR, which does not always cor- 
relate with the MIR. The 70 /im band of the Multiband 
Imaging Photometer for Spitzer (MIPS) instrument is 
not affected by PAH or silicate features for redshifts less 
than about 3, while the 160 /jm band is not affected at 
all. Moreover, these MIPS bands are closer to the FIR 
peak than the MIR data points used in previous studies. 
The longer wavelength MIPS bands should therefore be 
extremely useful for studying SMGs. The sensitivities 
and confusion limits at 70 /jm make this feasible, but dif- 
ficult, and hence only the deepest MIPS data are likely 
to lead to SMG detections. 

In this paper we present a study of the FIR properties 
and SEDs of submm galaxies using the deepest 70 and 
160 /im data available for the GOODS-N field. In par- 
ticular, we use the 70 and 160 /im data to check if the 
SEDs of distant SMGs are consistent with those of local 
ULIRGs. 

We assume a Hubble constant of 71 kms _1 Mpc _1 , to- 
gether with matter and cosmological constant density pa- 
rameters of rijvi = 0.27 and fl\ — 0.73 in this paper. We 
also use the notation SVo, Ssso, Si. 4 etc- throughout for 
the flux densities at 850 /im, 70 /im and 1.4 GHz. 

2. THE SUBMILLIMETER SAMPLE 

The SCUBA 'super-map' of GOODS-N contains 35 
robust 850 /im detections. Details of the data reduc- 
tion and source extraction can be found in Borys et 
al. (2003) and Pope et al. (2005). This submm image 
contains all publicly available SCUBA mapping data in 
this field taken up until 2004. Although 450 /xm data 



were also taken, they unfortunately have essentially no 
constraining power. Of the 35 SCUBA sources, 33 have 
likely radio and/or Spitzer counterparts using the Very 
Large Array (VLA), InfraRed Array Camera (IRAC) and 
MIPS 24 /im data in GOODS-N (see Pope et al. 2006). 
Thirty of these SCUBA galaxies are within the ultra-deep 
70 /im image area. SCUBA sources with low signal-to- 
noise ratios (SNRs) may have true flux densities which 
are boosted by a factor which depends on the source 
brightness and the local noise in the SCUBA map, and 
can be estimated if th e source counts are known (see 
e.g. ICoppin et al.l l2005). Therefore in this study we use 
the deboosted 850 /im flux densities given in Pope et 
al. (2006). 

3. SPITZER 70 AND 160 OBSERVATIONS 

The MIPS 70 /im observations were carried out during 
Cycle 1 of the General Observer program (Spitzer pro- 
gram ID 3325, Frayer et al. 2006). The inner 10' x 10' of 
GOODS-N was mapped to a depth of 10.6 ksec. The data 
were taken in the small-field photometry mode of MIPS 
with 10-second Data Collection Events (DCEs). Each 
Astronomical Observation Request (AOR) consisted of 
an 8-position cluster map, and the observations were 
completed with 12 AORs in total. In addition to our 
GO data, we used the MIPS G uaranteed Time O bservers 
(GTO) data (program ID 81. iDole et aT1l2004h . These 
GTO data have an integration time of 600 sec. 

The raw data were processed off-line using the Germa- 
nium Reprocessing Tools (GERT, version S1 3.1), following 
the a lgorithms derived by the MIPS team (|Gordon et al.l 
2005). Instrumental artifacts in the Basic Calibrated 
Data (BCDs) were removed using the filtering techniques 
adopted for the ext ragalactic First Look Survey (xFLS, 
iFraver et al.l [2006a') . The data were calibrated assum- 
ing an absolute flux calibration factor of 702MJysr _1 
per MIPS-70/im unit for stellar SEDs. The flux densities 
were multiplied by 1.09 to apply a galactic SED color 
correction. This assumes a power law SED of the form 
/„ oc v a and a = — 1, but the color corrections are sim- 
ilar for a — to —3. A more detailed discussion of the 
data reduction can be found in lFraver etail (|2006bD . 

The final image achieves a sensitivity of ^0.6 mJy rms 
and we have cataloged 101 sources (over 120arcmin 2 ) 
with SVo>2.3mJy (SNR>3<r). We catalog a region 
slightly larger than the area of deepest coverage to in- 
clude some relatively bright 70 /im sources and im- 
prove statistics. T he source counts are presented in 
IFraver et~a l. ( 2006b) and a full catalog will be presented 
in Huynh et al. (in preparation). The 70 /im image has a 
beam size of 18'.'5 FWHM, and in the presence of Gaus- 
sian noise the ler positional error of sources is of the order 
0.5 #fwhm/SNR, i.e. 3" for the faintest sources. 

The 160 /im observations of the GOODS-N region were 
taken as part of the MIPS GTO program in 2004. These 
data were taken in the scan mode of MIPS and we applied 
the standard filtering techniques to the 1 60 /im BCDs 
simila r to what was used for the xFLS (jFraver et alJ 
l2006a| ). The data were calibrated using a factor of 
44.7MJysr~ 1 per MIPS-160/tm unit. The 160 /im data 
have an effective integration time of 120 seconds and the 
160 /xm image reaches a sensitivity of ~ 15 mJy rms. A 
multiplication of 1.04 was applied to the 160 /im flux den- 
sities to color correct for galaxy SEDs. We also note that 
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the 160 fim light leak from 1 - 1.6 /im is not a problem for 
this data because there are no blue sources with mj > 5.5 
in the field. 

4. IDENTIFICATIONS AT 70 AND 160 /jm 

The negative K-correction at submillimeter wave- 
lengths means that SMGs are detectable in deep 
SCU BA images over a wide redshift range, 1 < z < 5 
(e.g. lBlain etai1l2002ft . However, at 70 and 160 fim, the 
K-correction does not compensate for the distance dim- 
ming, and the flux density of a galaxy with a given in- 
trinsi c luminosity drops s teeply as a function of redshift 
(e.g. lLagache et al.l [20051) . This is reflected in the high 
median redshifts of the SMGs (z ~ 2.0, Pope et al. 2006) 
compared to the 70 fim sources (z~0.5, Huynh et al. in 
preparation) . Hence the SMG and 70 /im samples are 
unlikely to have much overlap, and we only expect to 
detect the low redshift SMGs in deep 70 /im imaging. 

We examined the 70 fim image for counterparts to the 
SCUBA sources. To do this the 70 fim catalogue posi- 
tions were compared to IRAC positions of the SCUBA 
counterparts (Pope et al. 2006), searching within 10" of 
each submm counterpart. This search radius was chosen 
in order to take into account the typical positional un- 
certainties of low SNR 70 /im and SCUBA sources added 
in quadrature. 

This procedure uncovered two secure identifications of 
SMGs at 70 /im, GN26 and GN13 (following the nam- 
ing convention of Pope et al. 2005). These sources have 
70 /xm SNRs of 12 and 8, respectively, and hence 70 /im 
positional uncertainties of about l'/3 and l'/5 (including 
the Spitzer 1" pointing uncertainty). The positional off- 
set of the 70 fim source relative to the IRAC position 
is 0.6 and 0.3 arcsec for GN26 and GN13, respectively, 
which is well within the positional uncertainty. The prob- 
ability that one or more 70 fim sources lies randomly 
within a distance 8 of a SCUBA counterpart is: 

P = 1 - cxp(-7m6> 2 ) , (1) 

given a surface de nsity n of 70 fim sou rces (often called 
the P-statistic, e.g. lDownes et alll986| ). Hence the prob- 
ability of any 70 fim source lying within 1" of a SCUBA 
counterpart is 0.07%, using the 70 /im source density 
of 101 sources over 120 arcmin 2 , and hence the 70 fim 
matches to GN13 and GN26 are likely to be real. 

Two more SMGs (GN12 and GN32) have a nearby 
70 fim source at distances of 4.7"and 9.5", respectively. 
The 70 fim sources near GN12 and GN32 have IRAC and 
24 fim counterparts which do not match the identifica- 
tions for the SCUBA source (Pope et al. 2006). Although 
it is likely that some fraction of the 70 fim flux density is 
associated with the submm galaxy, it is difficult to deter- 
mine this fraction because of the other 24 fim sources in 
the vicinity. Using Equation (1), the random probability 
that 2 /30 SCUBA sources have a 70 fim source within a 
distance of 10" can be determined to be 28%. This is 
consistent with the 70 fim sources near GN12 and GN32 
being random matches. 

At 160 fim the only SMG detected is GN26. The 
beamsize at 160 fim is 40", so 70 fim and/or SCUBA 
sources which are close together may be blended into 
one 160 fim source. Examination of the 70 fim MIPS 
image suggests that the 160 fim flux density of GN26 
has some contribution from another 70 fim source (at 



z = 0.46) not associated with GN26. We therefore de- 
blended GN26 with a double Gaussian fit, fixing posi- 
tions to the IRAC counterparts of the two 70 fim sources 
which contribute to the 160 fim flux density. We find 
that GN26 has Si 60 = 110 ± 27mJy, which is 60% of the 
flux density of the 160 fim complex. The uncertainty 
at 160 fim is conservatively estimated to be 25%, taking 
into account absolute calibration, fitting errors, and de- 
blending issues. We also examined the SED of the source 
near GN26 to check the accuracy of the deblending. The 
5V0/S16O ratio of the second source at z = 0.46 is well fit 
by the quiescently st ar-forming galaxy SED templates of 
iDale fc Heloul (|2002fl . so the deblending at 160 fim seems 
reasonable. 

The multi-wavelength properties of the two SMGs 
which are detected at 70 fim (Figure 1 and Table 1) are 
described in detail in Pope et al. (2006). As expected, 
both sources (GN26 and GN13) are in the low redshift 
tail of the submm redshift distribution. Furthermore, we 
notice that these two sources are among the faintest at 
850 fim in the submm sample (both have Ss5o < 2.5 mJy) 
and therefore they are not typical of the full submm sam- 
ple presented in Pope et al. (2006) or indeed other sam- 
ples of SMGs. GN13 and GN26 have 70 fim flux densi- 
ties of 6.5 and 13.9 mJy, respectively, while the full 70 fim 
catalog has a median flux density of 5 mJy. About 80% 
of the 70 fim sources have spectroscopic redshifts, and 
the median redshift of these sources is z = 0.46 (Huynh 
et al. in preparation) . Thus the 70 fim counterpart to 
GN13 is typical of the full 70 fim sample, but GN26 is 
unusual in that it has a bright 70 fim counterpart, and is 
one of only seven 70 fim sources in our sample currently 
confirmed to be at z > 1 . 

5. STACKING ANALYSIS 

We performed a stacking analysis to derive an aver- 
age 70 fim flux density for the SMG population in the 
GOODS-N field. To begin with, we stacked the Spitzer 
data at the positions of all SMGs, including sources with 
70 fim matches or with coincident flux. For each SMG 
a square image 132" on a side (approximately 7 MIPS 
beams) was extracted. We rotated each image by 90° 
with respect to the previous one before co-adding to re- 
move any large-scale background effects. The median 
level of the individual extracted images was subtracted 
to remove any small scale offsets and yield better back- 
ground removal. Flux densities were determined using 
an aperture of 12" radius, and we applied an aperture 
correction of 2.0, which was calculated empirically from 
bright sources in the image. The aperture photometry 
was done after the background was subtracted from the 
stacked images, and the results were verified by making 
measurements with different size sky annuli. 

To estimate the expected scatter, offset stacked images 
were generated by randomly choosing a position in the 
70 fim image for each stacked source. Five hundred such 
random stacks were generated ,and the uncertainty in the 
stacked flux density is taken to be the standard deviation 
of these 500 measured values. 

The average 70 fim flux density ((6*70}) for all 30 SMGs 
is 2.00 ± 0.48 mJy, and hence the stacked signal is de- 
tected at over 4<r. The contribution from SMGs to the 
Extragalactic Background Light (EBL) at 70 fim can be 
estimated from this stacked signal by multiplying by the 
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appropriate source density. Assuming the SMG inte- 
grated source c ount of 2506 ± 406 deg~ 2 for Ssso > 2 mJy 
(jCoppin et ali r2006) , we estimate the contribution to the 
70 /im EBL from SMGs to be 0.016 ± 0.0 05 MJy sr" 1 
From a n extrapolation of the source counts iFraver et alJ 
(|2006bf) find that the total 70 /mi EBL is O.lSMJysr" 1 , 
so SMGs (5s5o > 2 mJy) make up about 9 ± 3% of the 
70 fim EBL. 

The 70 /mi stacked signal from all 30 SMGs is domi- 
nated by the 4 low redshift sources with coincident 70 /mi 
flux density. To determine the average properties of the 
SMGs not detected at 70 /mi, which is more representa- 
tive of the general SMG population, we also stacked the 
26 sources without coincident 70 /tm flux density. This 
time we stacked the residual 70 /mi image, which was ob- 
tained by removing all sources brighter than 3er at 70 /tm. 
Sources were removed by subtracting their fitted point 
spread functions from the image. The residual image 
was used to obtain a better signal to noise ratio in the 
stacked image and in the measured stacked flux density. 
For these 26 sources we find (6*70) = 0.70 ± 0.27 mJy, and 
so this stacked signal has an SNR of about 3. 

To test whether the majority of the stacked 70 /mi 
flux density is from the lower redshift SMGs, we also 
looked at low and high redshift sub-samples. The low 
redshift sub-sample consists of 12 SMGs out of the 26 
with z < 2, and the remaining 14 SMGs with z > 2 make 
up the high redshift sub-sample. For redshift complete- 
ness we included IRAC photometric redshifts from Pope 
et al. (2006) for 8/26 SMGs, and these all have z > 1.8 
(and estimated accuracy of Az < 0.4). The aperture flux 
density in the central region of the high redshift stack is 
0.22 ± 0.44 mJy, while for the low redshift stack it is more 
positive at 1.0±0.4mJy. This is consistent with the idea 
that the majority of the flux density from the full sample 
is coming from the lower redshift sources, although the 
measurements are too noisy to make a definitive state- 
ment. 

Similarly, we stacked the 160 /tm image at all SMG po- 
sitions, excluding the one detected source, GN26. There 
is no significant flux density in the stacked image and the 
3<7 upper limit is 13 mJy. A stack of the 12 low redshift 
SMGs gives a 3a upper limit at 160 /tm of 22mJy. 

To study the FIR properties of SMGs we could use 
the average flux densities of the full sample of SMGs. 
However, the stacked 70 /tm flux density from the sepa- 
rate high- and low- z sub-samples clearly show that most 
of the signal is coming from z < 2 sources (as expected 
from the K-correction). We therefore limit our analyses 
to the average properties of the low-z sub-sample in the 
following sections. 

6. INFRARED COLORS 

The average IR colors, Sjo/Ss^o and SVo/SW, are 
shown in Figure 2 as a function of redshift. GN26, at 
z~1.5 is consistent with a more active Dale & Helou 
(2002, hereafter DH02) model, with dust intensity index 
7=1.5, whereas GN13 has cooler colors, corresponding 
to 7 = 2.5. In the DH02 models 7 defi nes the amount o f 
dust as a function of heating intensity (|Dale et al.| [2001): 

dM(U) oc U-~<dU , 

where M(U) is the dust mass heated by a radiation field 
of intensity U. These 7 values span the range expected 



for infrared luminous galaxies - a low value around 1 im- 
plies a high contribution from photo-dissociation regions 
in an actively star-forming galaxy, while a significantly 
higher value represents the cirrus-dominated interstellar 
medium of a more quiescent or cooler galaxy. 

The infrared colors of the average low-z SMG are con- 
sistent with 7 ~ 2-2.5 DH02 SEDs. The average S m /S S r M 
ratio of SMGs from the stacking analysis indicates that 
the SMGs are relatively cool galaxies for their high IR 
luminosities, which is consistent with Pope et al. (2006). 
We also calculated upper limits to the S70 / Ssso ratio for 
SMGs not detected individually at 70 /tm (see Figure 2). 
The lower redshift SMGs are clearly inconsistent with 
lower values of 7. For the higher redshift sources we ex- 
pect much lower SVo/Ssso ratios, but with the current 
observational limits these sources can still be actively 
star-forming. 

For GN13, Figure 2 shows that the SJ0/S24 ratio is that 
of an actively star-forming galaxy, with 7~ 1-1.5, while 
the Syo/Ss5o ratio indicates it is a cooler galaxy. This 
may be due to broad silicate absorption falling into the 
24 /im band; the S70/S24 ratio can be strongly affected 
by PAH features and silicate absorption, which are not 
fully accounted for in the models. 

We find that the colors of GN26 are consistent with 
DH02 models with 7=1 and z~l-2 (Figure 3). The 
colors of GN13 place it at z—1-2 if a 7 = 1.5 model 
is assumed (Figure 3). However, GN13 is at z = 0.475 
and so this SMG has a warmer SVo/SW ratio than that 
suggested by its SVo/Ssso ratio, as mentioned earlier. 

The average colors of low-z SMGs are also plotted in 
Figure 3. We find that the average IR colors are best rep- 
resented with a DH02 model having 7 ~ 2-2.5 at z = 1-2. 
This is consistent with the median redshift of the low-z 
SMG sub-sample, and suggests the color-color plot can 
be used as a crude redshift indicator. 

7. SPECTRAL ENERGY DISTRIBUTION 

The FIR photometry at 70 and 160 /tm provides valu- 
able data points for constraining the FIR peak. Com- 
bined with the 850 /tm observations, the photometry 
spans both sides of the peak. Previous estimates of the 
SED of SMGs have relied on extrapolating the MIR or 
radio to fit the FIR peak - at the redshifts of SMGs, the 
MIR can be affected by complex emission and absorption 
features, so this method may not be reliable. 

We fit a variety of models to the data. These include 
four DH02 models with 7 = 1, 1.5, 2, and 2.5, as well as 
the Chary & Elbaz (2001, hereafter CE01) SEDs, which 
are templates derived from ISOCAM, IRAS and SCUBA 
data of nearby galaxies. The CE01 models have lumi- 
nosity dependent shapes, but since the local luminosity- 
temperature relation found for nearby galaxies may not 
hold for high-z SMGs (e.g. Pope et al. 2006), we fit CE01 
models allowing them to scale arbitrarily in luminosity. 

For GN13 and for the low-z SMG sub-sample, we con- 
strain the fit with the 70 and 850 /mi observations only, 
since they are not detected at 160 /tm (although the 
160 /tm data provide a useful upper limit). For GN26, 
the 70, 160 and 850 /tm observations are all used to con- 
strain the fit. We summarize the fitting results in Table 2, 
while Figures 4 and 5 show the best fit SEDs for GN13, 
GN26, and the average low-z SMG. Fitting CE01 SEDs 
without allowing the luminosity to vary freely results in 
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relatively poor fits, so we exclude these from the Ta- 
ble. This confirms the Pope et al. (2006) result that the 
local luminosity-temperature relationship does not hold 
for SMGs. For each best fit SED the total Lir (between 
8 and 1000 /im) was calculated and given in Table 2. Un- 
certainties in the luminosity were derived by scaling the 
best fit SED until the minimum \ 2 value exceeded the 
68% (±lc) confidence interval. 

The models provide good fits to the 70 and 850 /im data 
for GN13 and for the low-z average SMG, but they do not 
typically fit the 24 /im (observed) data point, although it 
has been shown that the 24 //m flux density can often 
be fit with additional extinction (Pope et al. 2006). The 
models do not provide a similarly good fit in the FIR for 
GN26; the 70 and 850 /im data points of GN26 are well 
fit, but the 160 /im measurement is underestimated by 
the models. Hence the luminosity for GN26 is probably 
higher than given by these models. Nevertheless, our 
derived IR luminosity for GN26 is 3 times greater than 
that estimated by Pope et al. (2006) from fitting SEDs to 
24 /im, 850 /mi and 1.4 GHz radio data. It is possible that 
there are further deblending issues at 160 /im for GN26 
(even although we have already divided the total 160 /im 
flux density between the two 70 /im sources in the area) . 
This demonstrates the power of 160 /im photometry in 
constraining the total infrared luminosity of galaxies, but 
shows that higher resolution is required to study such 
faint sources individually. 

The SMGs have high luminosities, but their FIR spec- 
tral shape is different from local ULIRGs of the same 
luminosity. We find that the average low-z SMG has a 
total IR luminosity of about 8.0 x 10 11 L©. This is a fac- 
tor ~ 2 less than the median SMG luminosity found by 
Pope et al. (2006) and Chapman et al. (2005) for their 
SMGs with z < 2. The reason that our calculated SMG 
luminosities are low compared to previous results is be- 
cause our best fit SEDs are cooler. The average SMG is 
best fit by a quiescent DH02 model with 7 ~ 2.5, or with 
CE01 SED templates of normal spiral galaxies scaled up 
by a factor ~ 300, with a rest-frame peak at about 150 /im 
(i.e. T~ 20 K). Local ULIRGs of the same luminosity as 
the SMGs are therefore not the best spectral templates 
for this sample. 

Several recent studies have relied on the MIR data (at 
24 /i m in particular) to derive luminosities and SEP fit s 
(e.g. iPerez-Gonzalez et all [20051 iLe Floc'h et~aLll2005[) . 
The 24 /im flux density was used in SED fitting of SMGs 
by Pope et al. (2006), who found that the typical SMG 
peaks at about 100 /an (corresponding to about 29 K). If 
the 24 /xm data point is included in the fitting of GN13 
along with the 70 and 850 /im data, we find the best fit 
DH02 SED peaks at shorter wavelength, corresponding 
to warmer dust temperatures, and the total IR lumi- 
nosity is decreased by a factor of about 2. There is no 
significant difference in the CE01 fits to GN13 with and 
without the 24 /im data point. For GN26 we find the 
24 /im data point makes no significant difference to the 
best fit DH02 SED, while the best fit CE01 model is 
slightly warmer and 3 times less luminous. The fit to the 
average SMG including the 24 /im data point is warmer 
than that with only the longer wavelength data, and the 
best fit DH02 and CE01 SEDs are about 2 times more 
luminous. At the median redshift of the average low-z 
SMG, (z) =1.4, PAH and silicate features fall into into 



the 24 /im band, and our fit here is driven by these fea- 
tures in the model SEDs. The DH02 and CE01 SEDs cer- 
tainly contain PAH features, but it is not clear whether 
SMGs at this redshift have strong or weak PAH features, 
if any. Therefore we would argue that the fit to the 850 
and 70 /im flux densities alone gives a more reliable re- 
sult for the average SMG total luminosity, at least for 
the moment, until we learn more about the MIR spectra 
of SMGs. 

8. DUST TEMPERATURES AND MASSES 

As a phenomenological alternative, we also adopt a 
modified blackbody SED model to fit the temperature 
of these SMGs. The SED is described by f v oc v B v , 
where B v (y, T) is the blackbody function for dust of tem- 
perature T, and /3 is the dust emissivity index. The 
MIR is approximated as a power-law of the form /„ oc 
v~ a and smoothly matches v^B v at longer wavelengths 
(Bl ain et al]|2003f h Although this simple phenomenolog- 
ical model cannot describe the full complex dust prop- 
erties of a galaxy, it can provide a good description of 
the general behavior of the SED. The range of parame- 
ters we consider is 15 K < T < 90 K and 1 < a < 4, which 
is representative of galaxies ranging from normal spirals 
to AGN. We set (3 =1.5 for our model fits, which is the 
yjdue found for dust in the galactic plane (|Masi et al.l 
Il995f). and a typical va lue for well-studied nearby galax- 
ies (jDunne et al.ll2000f ). 

We fit for T and a using the 70 /im and 850 /im data 
points for GN13 and the average low-z SMG, but also 
include the 160 /im detection for GN26. The results are 
summarized in Table 3. When fitting only two data 
points (and allowing the normalization to also be free) 
there is a strong degeneracy between a and T (which 
would be complete except for the boundaries of the pa- 
rameter ranges) - so the fit parameters must be inter- 
preted with caution. In the case of GN13, the full range 
of a is allowed by the 70 and 850 /xm data, with low val- 
ues of a corresponding to low values of T. Because of the 
additional 160 /im data point, the parameters are better 
constrained for GN26, with T~45K and a~3.5 being 
preferred. 

For the average SMG, the 70 and 850 /im flux den- 
sities alone cannot break the T and a degeneracy - a 
very low temperature of 15 K is allowed for a ~ 1.0, while 
T~33K for a~4.0. In a sample of 73 radio-detected 
SMGs the average 5*450 75*850 ra tio is measured to be 
5.0 ± 2.3 (|Chapman et al.l 12005). while 15 SMGs from 
this same sample have been detected with SHARC-II 
(Kovacs et al. 2006) and they have an average ^so/S'sso 
ratio of 4.0 ± 1.3. We find that models with a < 1.6 are 
inconsistent with these ratios. This implies that the al- 
lowed models for the average SMG have 21K<T<33K 
and 1.6 < a < 4.0, where the low values of T require low 
a. This shows that the low-z SMGs have relatively cool 
dust temperatures. 

The best fit dust temperatures of 21-33 K are con- 
sistent with those values previously derived for SMGs 
(e.g. Pope et al. 2006). Chapman et al. (2005) and Ko- 
vacs et al. (2006) suggested average temperatures close 
to the upper end of our acceptable range. This im- 
plies that the average low-z SMG in our sample has a 
relatively steep mid-IR SED, since our model fits with 
T ~ 30 K require a ~ 3. This suggests that the SMGs are 
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star-forming galaxies, because large a implies cool mid- 
IR colors, which are inconsistent with AGN-dominated 
sources. 

Assuming that the submm light is thermal emission 
from dust which is optically thin at A res t ~ 200 /im, with 
a single dust temperature T, the dust mass Md is given 
by: 

S&50 D T 



M d = r 7 S0U r , , (2) 

(1 + Z) Kd(^rost) B u (V Tcstl T) 

(e.g.lMcMahon 'et al.ll99J ). where Dl is the cosmological 
luminosity distance at redshift z and the dust absorption 
coefficient k is uncertain, even in the local Universe. We 
take a k value of 0.077 ± 0.030 m 2 kg _1 (jHughes et al.l 
1993), converting it to rest-frame frequency f res t with: 



Kd^rest) = 0.077 ( 



^rest 

350 GHz 



(3) 



Here we again assume that the dust emissivity index (3 
is fixed at 1.5. 

The range of allowable dust masses are calculated from 
the range of temperatures in Table 3. The dust mass cal- 
culated from Equation (2) is 1.0-1.6xl0 8 M Q and 2.2- 
2.5xl0 8 M o for GN13 and GN26, respectively, while 
the dust mass found for the average low-z SMG is 1.1- 
2.6xl0 9 M Q . This does not take into account the un- 
certainties in k and Sgso! the dust mass uncertainty is 
about 50% when these are added in quadrature. These 
dust masses are consistent, within uncertainties, with 
the molecular gas mass derived from CO observations 
(e.g. iFraver et all I1998L IFraver et all I1999L iNeri et all 
12003. lOreve et al.ll2005h . assuming a typical galactic gas 
mass to dust mass ratio M g /Md ~ 100 (e.g. iHildebrandl 
fl983h . 

9. FIR-RADIO CORRELATION 

Our sample can be used to test the FIR-radio corre- 
lation in submillimeter galaxie s. The FIR-radi o correla- 
tion is often expressed as (e.g. lYun et ai1l2001h : 



q = log 



FIR 



3.75 x 10 12 WnT 



log 



Si, 



Wm 2 Hz 



(4) 

where 'FIR' here refers to the flux between 40 and 120 
fim. The observed local value is q — 2.iA ± 0.3 (Yun, 
Reddy and Condon 2001). We use the best fit DH02 
models to derive the conversion from L40-120 to Ls-iooOj 
which is 2.0 and 1.6 for GN13 and GN26, respectively. 
Based on their Lir (Table [2]), we find q parameters of 
2.5±g;f and 2.4±° j for GN13 and GN26, respectively. 
These values are consistent with the local value of q, 
suggesting that these two sources follow the local FIR- 
radio correlation. 

10. CONTRIBUTION FROM ACTIVE GALACTIC NUCLEI 

As mentioned in Section 8, the inferred high values 
of a suggest that AGN do not dominate the bolometric 
luminosity in our sample of SMGs. To quantify the con- 
tribution of AGN to the infrared luminosity of SMGs, 
we adopt the simple modified blackbody approach as de- 
scribed in Section 8. For our AGN model we use f3 = 1.5, 
T = 90 K, and a = 1.1, as found for the xFLS AGN pop- 
ulation (Frayer et al. 2006). We subtract this AGN com- 
ponent from the observed 70 and 850 /im flux densities 



of the average SMG, and then repeat the fitting proce- 
dure with the DH02 and CE01 models, increasing the 
AGN component until the best fit % 2 value exceeds the 
previous minimum by la. This allows us to estimate the 
maximum AGN contribution to the FIR luminosity, with 
the results shown in Figure 5. 

It could be argued that the MIR waveband is the best 
discriminator of AGN, and therefor e we should be fo - 
cussing on the 24 /jm data point (e.g. lSajina et alJl2005f >. 
However, we are here calculating the percentage contri- 
bution of AGN to the Lir,, which is dominated by the 
FIR peak, and the contribution from AGN emission in 
the MIR is only a very small proportion of the total IR 
luminosity. 

For the average low- 2: SMG, an AGN which contributes 
up to 14% of the IR luminosity is allowed, using only the 
previous best fit CE01 SED model. If all CE01 models 
are used in the refitting, an AGN component of up to 
23% contribution is allowed. Together these fitting pro- 
cedures imply that the average (low-z sub-sample) SMG 
is dominated by a starburst. This is consistent with X- 
Ray studies which find that AGN contr ibute on average 
10% to the infrared luminosity of SMGs ([Alexander et al.l 
f2005h . 

For the 70 /im detected SMGs, GN13 and GN26, we 
find that 32% and 21%, respectively, of the total IR lu- 
minosity can be attributed to an AGN from the same 
analysis applied to the best fit templates. However, if 
the whole suite of CE01 SEDs are used, a larger pro- 
portion of the IR luminosity can be subtracted off the 
observed data points, with lower luminosity CE01 mod- 
els still fitting. This shows that with the current data a 
low luminosity starburst model from CE01 with an ad- 
ditional dominant AGN component is indistinguishable 
from a high luminosity CE01 ULIRG model, and that 
further photometry (or spectroscopy) is needed. 

Another method to estimate the AGN contribution is 
to make use of the 24 /im data. Here, we consider the ex- 
treme case where all the 24 /im flux density is due to an 
AGN, and determine its contribution to the total IR lu- 
minosity. For GN13 and GN26 such an AGN would con- 
tribute 11% and 5%, respectively, to the total IR lumi- 
nosity, assuming the T = 90K, a =1.1 AGN model used 
above, and the CE01 best fit SEDs. For the average low- z 
SMG, 21% of the total IR luminosity could be attributed 
to an AGN in this extreme case. This again supports the 
hypothesis that SMGs are dominated by star formation 
processes. 

11. CONCLUDING REMARKS 

We have presented 70 /*m properties of submillime- 
ter galaxies in the GOODS-N field. Out of 30 SMGs 
(with 5*850 > 1-7 mJy) in the overlap region of this field, 
2 are detected at relatively high significance in ultra-deep 
(~ 0.6 mJy rms) 70 /im imaging. Both of these detected 
SMGs lie at relatively low redshift. One SMG, GN26 at 
redshift 2=1.2, has infrared colors which indicate it is 
actively star forming. The second SMG, GN13 (z = 0.47) 
has infrared colors similar to normal spirals, but with 
a much higher luminosity. We confirm that these two 
SMGs detected at 70 /im follow the locally derived FIR- 
radio correlation. 

To determine the average properties of SMGs (most 
of which lie at z > 2 and are not detected individu- 
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ally at 70 /im) , we performed a stacking analysis and 
find that the average SMG has a 70 /im flux density of 
0.70 ± 0.27 mJy Most of the 70 /im flux density is com- 
ing from the lower redshift SMGs however. We analysed 
the average properties of twelve SMGs with z < 2, which 
have a stacked 70 /im flux density of 1.0 ± 0.4 mJy. From 
a stack of all 30 SMGs in the ultra-deep 70 /jm image, we 
find that the contribution of SMGs (with 5*850 > 2 mJy) 
to the total Extragalactic Background Light at 70 /im is 
9 ± 3%. 

The average low-z SMG ((z) =1.4) has cool infrared 
colors and an FIR SED that is best fit by a scaled up 
(about 300 x) normal spiral galaxy. We also found that 
an AGN contributes less than 23% of the total IR lumi- 
nosity in SMGs. 

We find that the average low-z SMG has an IR lu- 
minosity L 8 -iooo, of 8.0(±2.2) x 10 n L Q . The av- 
erage low-z SMG therefore has a star formation rate 
of 135 ± 35M Q yr~ 1 , using the relationship between 
SFR and IR lu minosity for starburst galaxies given by 
iKennicuttJ (| 1998ft . GN13 and GN26 have star- formation 
rates of 30 ± 10 and 800 ± 270M Q yr- 1 , respectively. 
The median IR luminosity of z < 2 SMGs in Chapman 
et al. (2005) is twice the value we find for our average 
low-z SMG, suggesting that theirs may have been an 
overestimate, suffering from lack of FIR data. 

The next generation submillimcte r bolometer, 
SCUBA-2, will come online in 2007 (H olland etall 



2006) and it is expected to yield a major improvement 
at 450 /xm, where it should reach a confusion limit of 
~ 1 mJy (5<r). This depth is well-matched to that of 
deep Spitzer 70 /xm imaging and should provide much 
more overlap than the 850 /jm selected sources. In 
addition, SCUBA-2 will produce much larger samples 
of fainter 850 /im sources, much like the two 70 (im 
detected submm sources in this study, and therefore 
70 (im data will be a valuable addition to understanding 
their SEDs. Also in the near future, the Herschel Space 
Observatory will produce confusion limited images 
across the wavelength range 75-500 /im. Herschel will 
therefore enable detailed study of the far-IR SEDs of 
large samples of SMGs, which we have shown here to be 
feasible for a sub-set of SMGs using Spitzer. 
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TABLE 1 

Summary of the properties of GN13, GN26 and the average low-z SMG. 



galaxy 


RA 


Dec 


s 24 


S 70 


s 160 


s 850 


s 1.4 






.12000 


.12000 


(f-Jy) 


(mjy) 


(mjy) 


(mjy) 


O-Jy) 




GN13 


12 36 49.72 


62 13 12.8 


371.0 ± 10.4 


6.5 ± 1.3 


<43 


1.9 ± 0.4 


45.4 ± 5.4 


0.475 


GN26 


12 36 34.51 


62 12 40.9 


446.0 ± 5.1 


13.9 ± 1.8 


110 ± 27 


2.2 ± 0.8 


194.3 ± 10.4 


1.219 


low-z SMG (12 soi 






258 ± 20 


1.0 ± 0.4 


<22 


4.0 ± 1.4 


<116 


1.4 



NOTE. — The 70 and 160 /^m flux densities include the absolute calibration error (of order 10%). The coordinates for GN13 and GN26 are IRAC positions (Pope et 
al. 2006). The redshift given for the low-z SMGs is the median for the sub-sample, and the Sqqq value given for the low-s sub-sample is the weighted average. 



TABLE 2 

Summary of the best fit SEDs and the calculated total Ltr. 



galaxy 


best fit DH02 


DH02 L m (L ) 


lum scaled CE01 Ltr (Lq) 


GN13 


7 = 2.5 


(2.5 ±0.3) x 10 11 


(1.8 ±0.2) X 10 11 


GN26 


7 = 1.5 


(4.5 ±0.7) x 10 12 


(5.0 ±0.7) x 10 12 


low-z SMG (12 sources) 


7 = 2.5 


(9.0 ± 2.5) x 10 11 


(8.0 ±2.2) x 10 11 



TABLE 3 

Allowable model parameters for a simple modified blackbody SED. 



galaxy 


T(K) 




a 


GN13 


34 


- 50 


1.0 


- 4.0 


GN26 


44 


- 48 


3.0 


- 4.0 


low-z SMG (12 sources) 


21 


- 33 


1.6 


- 4.0 



Note. — The range in T and a indicates the models which fit within the la limit of the 70, 160 and 850 fim observations. The dust 
cmissivity parameter, /3, is assumed to be 1.5. 
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Fig. 1— Greyscale HST ACS F606W images of GN13 (left panel) and GN26 (right panel). The contours show 70 /um data, plotted at 
3cr and 6a levels. The circle indicates the Pope et al. (2006) IRAC counterpart to the SCUBA source. The images are 40" X 40" in size. 
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Fig. 2. — Top Panel: SVo/SgSO c °l° rs °f the DH02 models as a function of redshift. From top to bottom, the solid lines are 7 = 1, 1.5, 
2, and 2.5. The dotted line marks an Arp220 SED, while the dashed line is Mrk231. The crosses mark the SMGs GN13 and GN26, which 
are detected at 70 (im. Other SMGs are shown as 3it upper limits. The filled circles denote the average ratios from the stacked analysis of 
the low-z sub-sample at (z) =1.4. Bottom Panel: S70/S24 colors of the DH02 models. Lines and symbols are as for the top panel. 
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Fig. 3. — S70/S24 versus SVo/Ssso colors compared with the DH02 models. For each model the squares mark redshift 0, 1, 2, 3, 4, and 
5, from right to left. The models from top to bottom are for 7=1, 1.5, 2, and 2.5. The crosses mark the SMGs GN13 and GN26, which 
are detected at 70 fim. The filled circle marks the average for the low-z SMG sub-sample. 
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restframe lambda (microns) 



Fig. 4. — SEDs of GN13 and GN26. The observed Spitzer and SCUBA data are plotted at rest-frame wavelengths. The best fit DH02 
model is plotted as a solid black line, while the dashed curve is the best fit CE01 luminosity dependent SED and the dotted line is the 
best fit CE01 SED with luminosity allowed to vary. The best fit Arp220 SED is also shown as the red solid line. The fits were constrained 
using only by the 70 and 850 (im data for GN13, and with the addition of the 160 fim data point for GN26. 
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Fig. 5. — The SED of the average low-z SMG. The observed Spitzcr and SCUBA data are plotted at rest-frame wavelengths and 
the curves are: best fit DH02 model (solid black line); best fit CE01 luminosity-dependent SED (dashed line); best fit CE01 SED with 
luminosity allowed to vary (dotted line); and best fit Arp220 SED (red solid line). The fits were constrained using the 70 and 850 fim data 
only, as described in the text. 
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Fig. 6.— The best fit SEDs for the low-z SMG sub-sample. The solid line is the DH02 model and the dashed line is the best fit CE01 
model (allowing the luminosity to float freely). The dotted line shows the AGN contribution allowed before all the CE01 models fail at the 
ler level. The symbols are as for Figure 5. 



